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Nonadiabatic tunneling in an ideal one dimensional semi-infinite periodic potential system is an- 
alyzed using our model. Using our simple analytically solvable model it is shown that ideal one 
dimensional semi-infinite periodic potential system can be thought as a model of moleular switch. 
The method is applicable to those systems where Green's function for the motion in the absence 
of coupling is known. The method has been applied to the problem where motion takes place on 
ON , parabolic potentials, for which analytical expression of Green's function can be found. 
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INTRODUCTION 



Quantum mechanical tunneling is an extremely important phenomenon in various fields of physics, chemistry and 
biology Most of the electronic transitions in molecular systems can be understood as non adiabatic transitions 
among electronically adiabatic eigen states at each fixed nuclear configuration. Nonadiabatic tunneling, is a tunneling 
accompanied by a nonadiabatic transition 0, S, 0] ■ The nonadiabatic tunneling is different from ordinary tunneling 
which involves single potential. It should be noted that these two cases become the same under strong diabatic 
Oh coupling limit. This is not the case in general, and the tunneling is affected by the upper adiabatic potential. So the 
w ■ nonadiabatic tunneling process can be controlled by the diabatic coupling, i.e., this can be switched on and off by 
switching the coupling. In this paper we are dealing with nonadiabatic tunneling in a one dimensional semi-infinite 
periodic potential system. Switching in one dimensional periodic potential systems have been gaining attention since 
the original works done by Kronig and Penney Complete transmission, which occurs at certain discrete energies 
in a semi-infinite periodic potential energy system, can be completely switched off by using a complete reflection due 
to an impurity potential created in the array. Positions of the complete reflections are determined by the shape of the 
i potential and the nature of coupling. So this complete switching can be realized by changing one of the potential in the 
t^. ' array to create an impurity. Similar idea was proposed by Carter Q, based on an array of periodic potential barriers. 
OO . The molecular switching model based on finite nonadiabatic tunneling potential system was originally proposed by 

(N ■ 

Nanbu et al. [7[ . In thise paper we use one dimensional semi-infinite periodic potential systems to show the possibility 
of molecular switching using our simple analytically solvable model. In the following we will first discuss our two state 
Cn ' model. 



OUR TWO STATE MODEL 



We consider two diabatic curves, crossing each other. There is a coupling between the two curves, which causes 
transitions from one curve to another. This transition would occur, in the simplest approximation near the vicinity 
of the crossing point. In particular, it will occur in a narrow range of x, given by 

V 1 (x)-V 2 {x)~V 12 (x c ) (1) 

where x denotes the nuclear coordinate and x c is the crossing point. V± and V 2 are determined by the shape of the 
diabatic curves and V\ 2 represent the coupling between them. Therefore it is interesting to analyse a model, where 
coupling is localized in space near x c . Thus we put 

V 12 (x) = K 6(x - x c ) (2) 

here Kq is a constant. This model has the advantage that it can be exactly solved 0,0,(13]. In the following we will 
first discuss our result of the two state curve crossing problem, which will be used in later section to analyze an ideal 
one dimensional semi-infinite periodic potential system. 
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TWO STATE CURVE CROSSING PROBLEM: EXACT ANALYTICAL RESULTS FOR THE DELTA 

FUNCTION COUPLING 

We start with a particle moving on any of the two diabatic curves and the problem is to calculate the probability that 
the particle will still be in that diabatic curve after a time t. We write the probability amplitude for the particler as 

Where rpi(x,t) and ip 2 (x,t) are the probability amplitude for the two states. *£(x,t) obey the time dependent 
Schrodinger equation (we put H = 1 ) 

d^(x.t) , s 

i—yj- = H*(x,t) (4) 



In the above equation H is defined by 



where Hi(x), H 2 {x) and V\ 2 {x) are 



H ~ 1 V 21 (x) H 2 (x) I (5) 



V 12 (x)=V 2 i(x) = V(x) (8) 

In the above V\{x) and V 2 (x) are determined by the shape of that diabatic curve. V(x) is a position dependent 
coupling function we shall eventually assume that V(x) is a Dirac delta function. In the following we will use the 
abstract operator notation of quantum mechanics. We found it convenient to define the half Fourier Transform ^>(E) 
of *(t) 

*(E) = / V{t)e lEt dt (9) 
Jo 

Half Fourier transformation of Eq. ((4]) gives 

^ 2 {E)) \ V E-H 2 ) [MO)) [ ' 

This may be written as 

= iG(E)$(0) (11) 
G(E) is defined by (E — H) G(E) = I. In the position representation, the above equation may be written as 

/oc 
G(x,x ;EjV(x ,E)dx (12) 
-oc 

The Green's function is defined by 

G{x,x ;E) = (x\(E - Hy^xo) (13) 
where \xq) denotes the position eigen ket with an eigenvalue Xq 

G(x,x 0] E)=(pf> X /> E J ] ^f,*o-,E)\ 
\G 21 (x,x ;E) G 22 {x,x ;E) J 
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Using the partition technique ll|we can write Gu(i, xq; E) as given below 

G u (x,x ;£) = (x\[E -Hi- V(E - H^V]- 1 ^) (15) 

The above equation is true for any general V. This equation considerably simplifies if V is a Dirac delta function at 
x c , which in the operator notation may be written as V = KqS = Kq \x c ) (x c \. Then 

G u (x,x ;E) = (x\[E-H 1 -K^G° 2 (x c ,x c ;E)S]- 1 \x ) (16) 

G 2 (x, xq; E) is defined by 

G° 2 (x 1 x ;E) = (x\(E-H 2 )- 1 \x ) (17) 

and corresponds to the propagation of the particle starting from xo, on the second surface in the absence of coupling. 
Now we will use the operator identity 

{E-Hx- K$G° 2 (x c , x c : E)S)- 1 = 

(E - Hi)- 1 + (E — H l )- 1 K->G° 2 (x c , x c - E)S[E -H x - K%G° 2 (x c , x c ; E)S}- 1 (18) 

Inserting the resolution of identity / = dy \y) (y\ in the second term of the above equation, we arrive at a 
Lippman-Schwinger type equation 



G n (x,x ;E) = G»{x,x ;E)+K£G 2 i {x c ,x c ;E) / dyG»{x,y;E)8{y - x c )Gu(y,x ; E) (19) 

J — OO 

G n (x, x ;E) = Gl(x, x ;E) + K 2 G?(ar, x c ; E)G° 2 (x Cl x c ; E)G 11 {x c , x ; E) (20) 
we now put x = x c in Eq. (|20p and solve it to find 

G n (x c , x ;E) = G°(x c , x ; E)[l - K$G%(x c , x c ; E)G° 2 (x c , x c ; E)}- 1 (21) 
This when substituted back into Eq. PU|) gives 

Lrn^afo,^ - i*x{X,X ,-U) + l-K*G a 1 {.x c ,x c ;E)G<H.x c ,x c ;E) (22) 

By using the same procedure one can derive analytical expression of G l2 (x, xq] E). 

r (~ ~ ■ F\ - K G°(x,x cl E)G° 2 (x c ,x ;E) (23) 
^12^; x 0, £>) — i-K%Gl(.x c: x c ;E)G%(x c ,x c ;E) 

Similarly one can derive the expression for G 22 (x, xq; E) and G 2l {x, xq] E). Once G(x, xq,E) is known, one can 
calculate wave function explicitely and using the wave function one can easily calculate transition probability from 
one diabatic potential to the other. In the following, we will apply our two state result to the case of one dimensional 
semi-infinite periodic potential system. 



ONE DIMENSIONAL SEMI-INFINITE PERIODIC POTENTIAL SYSTEM 

We imagine the problem of atom hopping reaction on a particular surface. Suppose active sites (which can adsorb 
atom) on the surface are arranged on a straight line. So we are thinking of one dimensional motion of an atom along 
the surface. The potential energy curve (after fixing the position of active sites) as a function of the position of atom 
will look like "one dimensional periodic potential system" . Now we consider the motion of active sites too. We take 

two neighboring active sites A and B, where an atom from A is transferred to B. The potential energy curve for A B 

* 

and A B crosses at a point (see Fig. [TJ. When this same atom transfers from B to C (the next neighboring active 

site), the potential energy curve for B C and B C crosses at another point. Here we assume that the motion of all 
active sites are coupled to the motion of it's nearest neighbouring active sites via Dirac delta functions - so that we 
can use our two state result to analyze this problem. So if an atom transports from one active site to another active 
sites through n active sites on the surface, there will be n + 1 potential energy curves and n crossing points (see Fig. 
[2]). Here we assume n to be very large (semi-infinite limit). In the following we will analyse this semi- infinite periodic 
potential system using our two state result. 



Co nf ig u ratio n a Coordinate 

FIG. 1: The potential energy curve for A B (left) and A B (right) crosses at a point. 
GREEN'S FUNCTION FOR ONE DIMENSIONAL SEMI-INFINITE PERIODIC POTENTIAL SYSTEMS 



Now first we assume that the all active sites are equivalent and the distance between any two neighboring active 
sites is the same, so the distance between any two neighboring crossing points will be the same. We have [n + 1) 
potential energy curves, we want to start our calculation with the n-th potential energy curve and let Gq (x, xq; E) 
denote the Green's function for the vibrational motion of the n-th active site, when it's vibration is uncoupled from 
the vibrations of all other active sites. Let G (ji ^ 1 \x,xq; E) denote the Green's function for the vibrational motion 
of first (n — 1) active sites, when only the vibrational motion of (n — l)-th active site is uncoupled from that of the 
n-th active site. Now using our two state theory (coupling is represented by Dirac delta function) we have derived 
an analytical expression of the Green's function for the vibrational motion of first n active sites, where the potential 
energy curve describing the vibrational motion of n-th active site is coupled to that of (n — l)-th active site at x n -\. 



G^(x,x :E) = G^(x,x ;E) + 



KjG^jx, x n - 1 ;E)G^- 1 Hx n - 1 ,x n - 1 ;E)G^(x n - 1 ,x ; E) 
1 - K^(x n - 1 ,x n - 1 ;E)G( n - 1 ){x n - 1 ,x n . 1 ;E) 



(24) 



Now we calculate the Green's function for the vibrational motion of the (n + 1) active sites, where the potential energy 
curve describing the vibrational motion of the (n + l)-th active site is coupled to that of the n-th active site at x n . If 
n is very large (semi-infinite) using the symmetry argument one can easily prove that 



G (B) (x n , x n ;E) = G*"- 1 ) {x n -i,x n - X ;E) 



So we get 



G { - n Xx,x ;E) = G%(x,x ;E) + 



K-0 Gq (x , Xn 



E) 



KlG™ + \x n ,x n ;E)G^)(x n ,x r , 



(25) 



(26) 
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FIG. 2: The potential energy curves as a function of vibrational motion of all active sites. 



If we put x = x n and xo = x n , in the above equation we get 

„ .T?\ — nn<„ ^ . i?\ i K oGo(xn,x n ~i- 1 E)G (n ^(x n ,x n ;E)Gn(x n ^i,x n ;E) 
(x n ,x n ,E)-G o{ x n ,x n ,E) + T^|Gp(^-^^ (27) 

The above equation can be simplified to 

-[KlG™ +1 (x ni x n ; E)G^(x ni x n ; E) - K^G^{x n ,x n - 1 ;E)GS(x n - 1 ,x n ; E) + (x n ,x n ;E) (28) 

K^ +1 (x n , x n ; E)G^(x n , x n ; Ef + G%(x n , x n ; E) = 

Now we solve the above equation for G^ n \x n , x n ;E). 

r(n)(„ „ ■ T?\ - [KlG^ + \x n ,x n :E)G^{x n ,x n -,E)-KlG^(x n ,x n . 1 -E)G'S{x n - 1 ,x n -E) + l] ft>Q . 
U {X n ,X n ,£.)- 2K%GZ + 1 (x n ,x n ;E) ^ 9 > 

, y / lK^Go +1 (x n ,x n ;E)G^(x n ,x n ;E)-K^G^(x n ,x n - V ,E)G^(x n - 1 ,x n :E) + l] 2 -4:K^Go +1 (x n ,x n ;E)G^(x n ,x n ^ 

2KlG^ +1 (x n ,x n \E) 

In the above expression, G^ (x n , x n \ E) is expressed in terms of Green's functions for uncoupled potential. For a 
continuum of vibrational eigen energies G^ n \x n , x n ; E) will have non-zero imaginary part for a range of E and that 
is possible only if the quantity A(E) as defined below, is negative. 

A = -4K^ +1 (x n ,x n ; E)G%(x n ,x n ; E) + [K^ +1 (x n ,x n ;E)G^(x n ,x n ;E) (30) 

-KlG n {x n ,x n - 1 ;E)G n {x n - U x n ;E) + l] 2 

So those enrgy ranges for which A(E) is negative, complete transmission is possible. In the following we will discuss 
the result of our numerical calculation using Eq. (30). 
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FIG. 3: The plot of "A" against E near the the energy of zoro-th vibrational state. 
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APPLICATION TO THE HARMONIC POTENTIAL CASE 

It is possible to find Gq (x, Xq]E) only for few cases and harmonic potential is one of them. For example, if 



(31) 



then the analytical expression of the corresponding Green's function is known [12j. We calculate A for harmonic 
potentials [Eq. (30)]. We use harmonic potential of frequency uj = 500 cm' 1 and mass m = 35.4. The distance 
between the equilibrium positions of any two neighboring harmonic potential we take, is 0.1 A The distance between 
the crossing points of any two neighboring harmonic potential we take, is 0.05A. K we take is 1.58 x 10~ 7 erg. A 
The result of our calculation (see Fig. [3]) shows that "A" is negative over a range of energy near zero point vibration 
(t^), this means that G^ (x n , x n \ E) has non-zero imaginary part over the same range of energy, so the system has 

a continuum of vibrational eigen energies in that range of energy. Same is true for other vibrations (^f 2 -, ^j 2 - ) 

also. So complete transmission is possible near those energy values ^p-, ^f 1 ). These complete transmissions 

can be completely switched off by using a complete reflection due to an impurity potential created in the array. 



CONCLUSIONS: 



Nonadiabatic tunneling in an ideal one dimensional semi-infinite periodic potential system is analyzed using our model 
Hamiltonian. It was shown that complete transmission is possible for certain range of energies and complete reflection 
is possible by adding impurity potential in the array. The same model can be applied to multi-state electron transfer 
model, where the atomic coordinate will be replaced by a global coordinate representing the nuclear polarization of 
the system. Our method is applicable to those systems, where Green's function for the motion in the absence of 
coupling is known. The method has been applied to the problem where motion takes place on parabolic potentials, 
for which the Green's function can be found. The same procedure is also applicable to the case where S is a non-local 
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operator and may be represented by S = \f) Kq (g\ , f and g are arbitrary acceptable functions. Choosing both 
of them to be Gaussian should be an improvement over the delta function coupling model. S may even be a linear 
combination of such operators. 
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